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A USER’S GUIDE FOB V174 
A PROGRAM USING A FINITE DIFFERENCE METHOD 
TO ANALYZE TRANSONIC FLOW OVER OSCILLATING WINGS 

T. D. Butler, W. H. Weatheriil, 

J. D. Sebastian, and F. E. Ehlers 
Boeing Commercial Airplane Company 

1.0 SUMMARY 

This document describes the design and usage of a pilot program for calculating 
pressure distributions over harmonically oscillating wings in transonic flow. The 
procedure is based on separating the velocity potential into steady and unsteady parts 
and linearizing the resulting differential equation for the unsteady flow by assuming 
small disturbances. The steady velocity potential distribution, which must be obtained 
from some other program, is required for input. The differential equation for the 
unsteady flow is linear, complex in form, with spatial varying coefficients that are 
dependent on the steady velocity potential distribution. Time is not a variable in the 
program since sinusoidal motion is assumed. The numerical solution is obtained 
through a finite difference formulation and a line relaxation solution procedure. 

This program may be used for wings with arbitrarily, swept leading- and trailing-edges. 
The leading edge may be curved. However, the trailing edge must be straight. The 
program uses a rectangular array of mesh points. 


2.0 INTRODUCTION 

This document describes the design and usage of the FORTRAN IV digital computer 
program V174. This pilot program was written as part of a research effort to develop a 
method of computing the transonic perturbation flow about a harmonically oscillating 
three-dimensional wing. 

Based on the finite difference procedures of references 1 through 3, the program 
computes the solution of the unsteady velocity potential and the resulting unsteady 
pressure distributions. It requires as input the potential distribution from a steady state 
transonic small perturbation program (ref. 4). Conservative differencing is used for 
subsonic points and nonconservative differencing for supersonic points and across 
shocks. Figure 1 shows the geometry used for the program. 


The program is set up to calculate the deflections and slopes for a control surface mode, 
a pitch mode, or a flapping mode. The control surface, which is located at the trailing 
edge, may be either partial or full span. Although a modal amplitude theta (0) is input, 
the printed pressure coefficients are for unit amplitude motion. Other motions may be 
treated by modifying the appropriate boundary condition subroutine. 






Included in this document are: 
e Description of equations used in the program 
e Description of user I/O and scratch file formats 
e List of program limitations 
« Description of computer program usage 

a Data stacking procedures 

I 

a Description of program output, normal and diagnostic 

9 Sample problem input/output ■ - 

a Description of program structure and routines used 
Features of this program include options to: 
e Analyze flat plates or wings with thickness 
a Use row or column line relaxation 

a Begin the iterative solution from a previous solution of the same dimensions 
a Analyze wings with swept leading and trailing edges 

a Take advantage of autosymmetry in the case of a wing at zero angle of attack and 
with a symmetrical thickness distribution to substantially reduce the number of 
nodes 

a Use over- and under-relaxation factors to reduce the number of iterations to 
convergence 

a Determine iteration print intervals 

a Determine far-field updating intervals 
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3.0 SYMBOLS AND ABBREVIATIONS 


a,b 

coefficients for y,z differences corresponding to second deriatives, 
with appropriate subscripts (eq. (3)) 

b 

root semichord of wing 

c,d 

coefficients for x difference corresponding to second derivative 
(eq. (3)) 

Cki.Ck2iCk3>Ck4 

equation (28) 

Csl>*'s2>dsl)ds2 

equation (13) 

ci,di,C2,d2 

coefficients for second-order accurate difference corresponding to 
first derivative (eq. (3)) 

E 

coefficients in difference equations with appropriate subscripts 

Fij 

equation (6) 

f(x,y,t) 

instantaneous wing shape defined by zq = 8f(x,y,t) 

fo 

undisturbed wing or airfoil shape 

h 

unsteady contribution to wing or airfoil shape 

h 

Xm+l 

i 

VT 

i.j.k 

x,y,z subscripts for points in the mesh 

M 

freestream Mach number 

U 

freestream velocity 

Xa 

coordinate of control surface hinge line 

x,y,z 

scaled coordinates (xq, /iyo> /^Zq) for the three-dimensional problem 

x,y,z 

coordinates defined in equation (18) 

Xjj.Xt 

coordinates of leading and trailing edges 

Xo.yo.Zo 

physical, coor,dinates, made dimensionless with the root semichord 

Xi',yr,zi' 

variables of integration 
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coordinate of wingtip 


yt 

/3 

y 

ACp 

AX] 

Axg 

A^i 

S 

8l 

§2 

e 

^1 

<t> 

<Po or Phio 
tpi or Phij 
<Pl 

w 

X 

0 ) 


Vl-M^ . 

ratio of specific heats for air • 
jump in pressure coefficient 

Xi - Xi.i 
Xi_i - Xi_2 

jump in at plane of wing or vortex wake 

thickness ratio or measure of camber and angle of attack 

X2 -xi 

^Hnax ■ 

(8/M)^^^ 
wM/(l - M^) • 

scale factor on and Zq, (l = 8^^^ 
unsealed perturbation velocity potential 

steady scaled perturbation velocity potential 
unsteady scaled perturbation velocity potential 

wake integral defined in equation (21) 
acceleration or pressure potential 

fundamental source solution of integral equation for evaluation of 
far-field boundary conditions 

angular reduced frequency (semichord times frequency in radians 
per second divided by the freestream velocity, wbAJ) 
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4.0 DISCUSSION 


The finite difference equations for three-dimensional unsteady transonic flow are given 
by Ehlers in reference 1. At interior points where the steady flow is subsonic, we have 

and at points where the steady flow is supersonic, 

"^k ijk-l ■ (“''i * '’S'j ^ ijk ■ (2) 

*’% ‘'’■ijk+1 ° ^ ’’■i-ljk ■ i-2jk ' '^'S-lk ' '’''j '^‘ij+lk 


where: 


1 


^k (^k+l ■ ^k-l^^^k ■ ^k-P 


^yj (yj+i -yj-i)(yj-yj-i) 

1 


’^k (Zk+I ■ ^k-1^ ^^k+1 ‘^k^ (yj+1 ■ yj-lHyj+i -Yj) 


E. = 


E. = 


E, = 


E. = 


"i^i+yzjk 

^i+'/zjk + 

^i-1 

di-lUi-3/2jk-’"d2i/^ 
1 


(Xi+i - Xj.p (XJ+.2 - Xj) ^4 ^1-1^ 


i+ljk ^°ijV 


^i+Vzjk" K - (7 + l)(‘Po;^i;i, ' 


«i-y.jk= K - (T + l)(v^„yj^ - 


(3) 



--1 

5^1-2 ■ 


Xj-Xj.l 




-di. 1 
■^1 *1-1 


Equations (1) and (2) are equations (24) and (27) from reference 1. 

The boundary conditions on the upper and lower wing surfaces lead to the following 
equations for subsonic flow at finite difference points immediately below the wing, 
k = k, 


m 


'■ (“yj * S ■ “"k. ■ ■"= 


^2.1 •• 


+ a, + E, + E, 


ijkj„ 


(4) 


■El i+ijk.- '^2 'f'l i-ntf __ '^y,- *^1 *^1 ii+lk_ 'El El 


(L)- 


m 


i-ijk^ ^yj^iij-ik^ -yj-iij-Mk^ - 


and points immediately above the wing, k = km + 1 

'’■i+ljk^+1 ■^"^li-lik^+l '“yj *’>ij-lk„+l ■'■ij+lk^+1 ^ 

where 

Fif 1 = (Xj. yj) + to f;i-l (Xj, yj) 


f.™ <ki. yj) + i" f.‘“* (kj, yj) 


The (L) and (U) refer to upper and lower wing surfaces, respectively. The equations 
similar to equations (4) and (5) at supersonic points can be written down analogously. 


The total harmonic deflection of the wing is written as 

z„ - 5f(x,y,t) ^ 6 {f„(x,y) + (x,y) } 


(7) 


The steady velocity potential (^iq) is calculated from the steady deflection shape 
(fo),while the unsteady potential {<p{) is calculated from the harmonic mode shape 
fi(x,y). 


Over the wake, the condition that the trailing vortex sheet supports no pressure, 


dAtPj 

3x 


+ icoAv2j = 0 


(8) 


IS 

OF PCX)R QUALITY 



results in a term being added to the right-hand side of equations (1) and (2). For finite 
difference points just below the wing plane (k = km), the additional term is 


*^m y 


(9) 


and for points just above the wing plane (k = km -f- 1), the term is 


where 


-a„ A<A .. 

^kn,+ l 


A<p, = Av? e ■ ^i,+l^ 
Mj + lj ' 


(10) 


( 11 ) 


and is the jump in velocity potential at the first point aft of the wing trailing 

edge at station j determined so as to satisfy the Kutta condition on the trailing edge. 
The addition of equations (9) and (10) implicitly satisfies the conditions. The normal 
velocity is continuous across the wake. 

The finite difference equation for the jump in across the wing to the second order in 
mesh size is 

A.. - 


( 12 ) 


where 


d.i - 


s, = 


1 

1 


4sj(Sj + 1) 

""sa 4s, (s, + 1) 


h(2Sj + 1) 
4(Sj + 1) 

h(2s- + 1) 
1 ^ 

s 2 4 ( 82 + 1) 

(13) 

,+2 -^k^+l)/'’ 




h = Zi, -Zy. 

‘ m 


(14) 
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Two integral relations are used to satisfy the far-fleld boundary conditions on the outer 
boundaries of finite difference mesh. The first for the velocity potential is 





+vt 

J ./ . t""' 

-Yt xg(y,') 




iojXtCy/) , , f 

e )dyj J 

Xt(yi') 


-icoxi 

e 


dx,' 



(7+l)<Po .<^1 , 

Xj Xi 




(15) 


and the second for the pressure function is 

P(Xj,yi,Zj) = 'f I + 


(16) 


_L 

4 n 


1 

'W 


+Yt xt(y;) 

■ j [ 

-Yt X£(y/) 

+yt 


-yt 


-XAv^j ,1dx;dy/ 
Z- -J 


^ (xt(y/)-x,,yi -yi',Zj)dy/ 


^4jr| Xx,' ’'^•x.'x;)"''' 


(17) 


defining 


= Xj -x/, y = Yi -y,', z = z, - z/ and R 

giXj(Mx-R) 


I 2 3 ^ 
=V "iff + y z 


i/^Cx^y^z) = 


R 




'/'x; = -^(i/R+i\)]'^ 


(18) 
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( 18 ) 


X = li/v + 

.1 

X = [ iX,M -^(1/R + iXj ) + iw] i// 

X^' = -^|[^ - iXj M - icoj (iXj + 1/R) - x/ x/r| i/' 

Equations (15) and (17) have been simplified for the program. First, as noted in the 
two-dimensional derivation 


9Ai^j 

9zj 



= 0 


(19) 


and thus the second integral in the first term of both equations (15) and (17) is zero. 
Second, the third term, which is the volume integral and has not been of significance in 
the two-dimensional problem, has been dropped. Third, since we are interested in the 
far-field, we approximate Xj - xj' and yj - y^' with yj so that the terms of and x may 
be moved outside the integral sign. The evaluation of wake integral in equation (15) is 
discussed in detail in the next section. The equation for the velocity potential on the 
far-field (eq. (15)) for Xg = 1.0 is 


</>i(x,,yi,z,) - 



dx.'dy,' + 'P. 
Ill iw 




where 


'Pi 


w 


+yt 


( 21 ) 


-yt 


with the trailing edge defined as a straight line by the functions f(yi') - a| yi‘ |, the 
function is given as 

^ - f , / U] 






1 


R 


0 


[y l+ui^ \ /3yf2(yi') + Ro^ RO 


+ ■ 


i^ 

M 


N 


■p yi + ui2+ ^ WjF(Tj) 

i = 1 


( 22 ) 
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Ro=y(yi + 


The pressure function equation (9) becomes 

+yt 

P(Xj,yi,Zi) = j J A<PjdXjdyi 

-yt 4(y/) 


;yt 

1 e Wy/) A^. (y/)dy; 

4tt J ‘ ‘ 


Equation (20) is used to evaluate the velocity potential along the line resulting from the 
intersection of the Y-Z plane through the trailing edge of the wing and the X-Y and X-Z 
planes bounding the finite difference volume. Equation (16) is then integrated by the 
trapezoidal rule to determine values ahead of this line and behind this line on the 
upper, lower, and side boundaries. For example, on the lower boundary when k = 1, for 
points ahead of the trailing edge (x<1.0 or i<ii), 

yP(xi,y;,Zj)ei‘^(^i'>^i-l) + P(Xi.i,y.,z^^^ 

M-ljl ‘ijl 

and the equation for points downstream of the trailing edge (ii<i) is 

+ [p(xj,yj,Zj) +P(xj.j yj,Zj) e 

The application of equation (23) to the upstream and downstream boundaries results in 
the following equations: on the upstream boundary, 

(26) 
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and for the downstream boundary, 


jk '^1 j ] jji ‘^k4 ijnaxJ^ 

'max ■’ 


Pljk - P(Xi,yj,zi,) 


Pi ik " P^^i 

‘maYJ''- ‘max’ J 


1 + icoS j /2 
and Ckl 1 - ico5j/y 

Ck-2 = 5i/0 -IOJ8J2) 


5i = X2 -Xi 


0,3 = (l-ico52/2)/n +ico52/2) 63 = - 1 

Ck4 - 52/(l+ico52-/2) 

Equations (25) and (26) may be used to substitute for and in equations (1) 

and (2). 
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5.0 COMPUTER PROGRAM USAGE 

5.1 MACHINE REQUIREMENTS 
V174 executes on a CDC 6600 or similarly compatible computer. 

5.2 OPERATING SYSTEM 

V174 was'designed for a KRONOS 2.1 or NOS 1.1 operating system. 

5.3 STORAGE ALLOCATION 

This program executes in 132000 octal words of computer memory. 

5.4 TIMING 

Timing is hardware and operating system dependent. 

The following example was run on a CDC 6600 with KRONOS 2.1 operating system. 
The program loads in about 4.5 CP (central processor) seconds. Using an XYZ mesh of 
44 by 16 by 26 (18 304 nodes), the program requires about 5 CP seconds per iteration. It 
needs a like amount of time for each far-field update; 50 iterations of this example with 
a far-field update every 10 would require: 

CP seconds = 4.5 + (50)(5) + (50/10 + 1) (5) 

or 284.5 decimal seconds. Rounding up the user would place 300 decimal seconds on his 
job card. 


5.5 FILE I/O 


The program card of V 174 is as follows: 

Program V174 (INPUT= 1002,OUTPUT=1002,TAPE7, DEBUG, TAPE5=INPUT, 
TAPE6=OUTPUT,TAPE1=1002,TAPE10=1002, 
TAPE11=1002,TAPE12=1002,TAPE13=1002) 

As noted previously, the buffer sizes of most of these files have been reduced from the 
default values to save memory space. 

5.5.1 FILE UTILIZATION 

A considerable saving of core membry, and consequently cost, was achieved by 
modifying the program. This modification puts only 3 planes of the Phii and PhiQ 
matrices in core at a time. These planes are chordwise X-Z planes and are denoted by 
the J indice (identifying their Y coordinate). Thus, on any one binary file there will be 
JMAX planes of data. 


ORIGINAL TAGE IS 
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In core the program identifies the necessary 3 planes as JR (corresponding to J+ 1-see 
theory), JC, (corresponding to J), and JL (corresponding to J-1). There are 3 indices 
available in the incore matrix data space, but each of these will be reused as often as 
necessary moving JR, JC, and JL. For instance, if the program is calculating the X-Z 
chordwise plane corresponding to J indice 6 this plane may have been read into the 
incore storage area plane numbered 3 and would designate that incore plane as JC. 
Calculation plane corresponding to J = 5 would be in JL which would be stored in the 
incore data storage plane 1. The next calculation plane JR (corresponding to J=7) would 
be read into the available incore plane storage area 2. 

After calculation plane J = 6 is recalculated, and written to binary file, this plane is 
now designated JL by the program. The former JR now becomes JC (incore plane area 1) 
and a plane corresponding to J = 8 is read in and placed in the incore storage area now 
designated JR (incore storage area 3). 

Thus, for calculation of plane J = 6 the incore storage should be as follows; 

Calculation Incore plane 

plane storage area 

Plane JC 6 3 

Plane JR 5 1 

Plane JL 7 2 

For calculation of plane J = 7: 

Plane JC 7 1 

Plane JR 6 ^ 

Plane JL 8 3 

Please note that the description applies to the usage of the Phio matrix as well, with the 
exception that it is not recalculated with each iteration. 

The file named TAPEl contains one steady state velocity potential distribution (Phio) 
from a separate program. This matrix must be present for a wing with thickness (see 
ref. 3 and input variable list MSTST). 

The file named TAPE13 contains the modified Phio matrix from file TAPEl. l^is 
matrix is written on the file in the standard form for this program. The initialization 
overlay recalculates this matrix. During iterations the Phio matrix on this file is read 
as previously described. 
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The file named TAPEIO is the file containing the starting data, if any (see input list 
variable, INC), for the Phij matrix. After a. successful execution of the program it will 
contain the new Phij matrix in place of the old. Its format is standard for this program 
,(sec. 3.5.2). 

The files named TAPEll and TAPE12 are interchanged in function for each iteration of 
the program. They will contain the scratch Phi^ matrix in the standard format 
discussed in section 3.5.2. Initially, one of these will contain the starting Phij matrix. 
Either may, be from TAPEIO, a previous run, or a Phij matrix of zeros. The other will 
then have written on it the Phi^ calculations from this iteration. For the ngxt iteration 
these files will be rewound and functions interchanged with the former supply file 
becoming the receiving file for this iteration. 

5.5.2 FILE FORMATS 

Binary files 

The files TAPEl, TAPEIO, TAPEll, TAPE12, and TAPE13 are binary I/O files and all 
use the same basic structure of format. TAPEl and TAPE13 are real data (Phio matrix), 
TAPEIO, TAPEll, and TAPE 12 are complex with a real- and imaginary -parts written 
for each element of the matrix (Phij). 

The structure of these binary files are as follows: 

1 (J indice of the first plane) 

The l'^^ element at K = 1 


The IMAX^^ element at K = 1 
The I^*^ element at K = 2 


The IMAX^^ element at K = 2 


The IMAX^*' element at K = KMAX 


2 (J indice of the second plane) 

• same structure as the above plane 



JMAX (indice of last plane) 
element at K = 1 


IMAX*'*’ element at K = KM AX 

As mentioned previously the complex binary files will contain a real part and an 
imaginary part for each element of the Phi^ matrix. 

BCD Files 

The BCD files INPUT, OUTPUT, and TAPE? follow the standard FORTRAN and 
system formats for that type of file. 

BCD files are those which deal with character printing or reading. V174 has three 
of these; INPUT, OUTPUT, and TAPE?. 

INPUT, also called TAPES, is the file which contains cards or card images. Program 
card (or card image) input is fully described in section 3.7. 

OUTPUT, also called TAPES, is the file on which the program places the primary 
printed information. (See sec. 3.8 and app. A.) 

i , t ‘ , 

TAPE? is also a, print file. The user may disregard it unless he is executing the program 
on a terminal where the primary OUTPUT print file is to be printed later. It will print 
a summary during execution telling the user how convergence is proceeding. 

Usage of TAPE? is LGO,IN,OUT,OUTPUT. (See sec. 3.6.) 

TAPE? was primarily used in development of the program. Terminal usage of the 
program should be limited as terminal execution is usually very expensive. 

5.6 CONTROL CARDS 

The following control cards can be used to load and execute TEV174 from tape; 
JOBN,T500,CM132000,P02. 

ACCOUNT,ACCTNO,PASWD. URNAME/PH/M.S./ORG 
REQUEST, TAPE,VSN=66XYYY,F=I,LB=KL. 

COPYBF,TAPE,LGO. 

RETURN,TAPE. 



GET.TAPE1=URPHI0. 

GET.TAPE 10= URPHIl . 

MAP, FULL. 

LGO. 

End of Record 
Input Cards 

The preceding will execute using existing Phigand Phii distribution files. If either is not 
available, then adjustments are made in the input cards and the respective GET card is 
not used. 


5.7 PROGRAM INPUT 


5.7.1 GENERAL REMARKS 

The input to TEV174 is of two forms, disk file/tape input (binary input) and card input 
(BCD input). Disk file/tape input may consist of input point locations and the Phio 
and/or Phij distribution. An input Phij distribution is indicated if a previously 
calculated Phii matrix is to be used to start the iteration process. If the user is starting 
from scratch, there is no Phij input and the initial Phij distribution is ail zeros. Phio 
the steady-state distribution from another program. If a flat plate solution is. sought, it 
would not be input. 

A description of disk file/tape formats is given in section 5.5.2 and a listing of the 
input for a sample problem is presented in appendix A. 

The card input consists of field dependent input and namelist free field input. The field 
dependent input is defined in the format column of table 1 as a specific field (i.e., F10.2, 
AlO, 15). The namelist data will be represented in the same column by the namelist 

name "PARAM.” 

Some of the features of namelist input are: 

1. Card field consists of columns 2 through 80. 

2. List consists of a $ list-name in column 2 followed by a series of specifications 
continued on as many cards as required and terminated by a $. 

3. Specifications are of the form: 

a. Vname = Value 

b. Where Vname is an array, Vname = Valuel, Value2, . . . ,Valuen 

Where Vname is one of the variable names for the list, value is the associated 
value(s). Value may be an integer, a floating point number in normal or 
exponential form, or in the case of a logical variable (specifically the options) of the 
form. 
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.T. or .True, indicating true or on 
.F. or .False, indicating false or off 


4. Specifications must be separated by commas. There is no comma between the last 
specification and terminating $. 

5. Embedded blanks are allowed except within the $ list-name, variable name, or 
value. At least one blank must separate the $ list-name and the first specification. 

6. The order of appearance of variables on the card(s) is not important; the spelling is. 

7. Any or ail of the variables may be left out of the list, e.g., $ list-name . . $ is 
legitimate. This assumes, of course, that there is a legal default value associated 
with the variable(s) not included in the list. 

5.7.2 LIMITATIONS 

The following are size limitations within the program. 

9 3^IMAX^55 X nodes parallel to upstream downstream flow in the XYZ 

mesh. 

9 3^JMAX=s 33 Y nodes spanwise in the XYZ mesh. 

9 3«KMAX^26 Z nodes vertically in the XYZ mesh. 

Note; The Phii and Phio distributions must also correspond to the limitations on 
the XYZ mesh. 

9 l=sNDWING^20 The number of Y values in the Y array less than YT. 

Corresponds to number of values in XLE, XA, and XT. 

5.7.3 DATA STACKING 

Note: All coordinates are entered as scaled coordinates, 

X = Xp/b 

Y = /i * Yp/b 

Z == fi* Zp/b 

where b is the root semichord and the subscript p means physical coordinate. 
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Table 1. - User Input Variables 


Card number 
1 

2 to N 


Variable name 

Format 

Description 

Title 

8A10 

80-character title for this problem 

FSMACH 

PARAM 

Mach number 

DELTA 

PARAM 

Thickness ratio 

THETA 

PARAM 

Maximum angle of attack in degrees 

OMEGA 

PARAM 

Angular reduced frequency 

GAMMA 

PARAM 

Ratio of specific heats for flow medium 

AL 

PARAM 

Axis of rotation for flapping mode 
dimensionless physical coordinate 

IMAX 

PARAM 

Maximum X node count in users mesh 

JMAX 

PARAM 

Maximum Y node count in users mesh 

KMAX 

PARAM 

Maximum Z node count in users mesh 

IS 

PARAM 

Starting X node limit for volume integral 
for the wing; currently not used 

IE 

PARAM 

Ending X node limit for volume integral for 

r 

the wing; currently not used 

KS 

PARAM 

Starting Y node limit for volume integral 
for the wing; currently not used 

KE 

PARAM 

Ending Y node limit for volume integral for 
the wing; currently not used 

NMAX 

PARAM 

Maximum number of iterations to be 
allowed without convergence or divergence 

NA 

PARAM 

Far-field update cycle control; updates the 
far field each NA iteration 

ERROR 

PARAM 

Error difference. - When the maximum 
difference between Phii distributions of 
consecutive iterations is less than ERROR, 
the program stops iterating 

NP 

PARAM 

Prints pressure distribution every NP 
interations 

INC 

PARAM 

Restart variable. - If INC=0, start with the 
Phil distribution on TAPEIO; if INC=0, 
start with a Phii distribution of zeros. 

ORF 

PARAM 

Overrelaxation factor used for subsonic 
nodes to accelerate convergence 1.^0RF=s2. 

URF 

PARAM 

Underrelaxation factor used for supersonic 


nodes to accelerate convergence 0.<URF^1. 


OF FOOB QUAIXTY 
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Table 1.- (Continued) 


Card number Variable name Format 
2 to N (cont) MSTST PARAM 


ISWEEP PARAM 


ILAX PARAM 


CONPXT 

PARAM 

CONE6 

PARAM 

YS 

PARAM 

YSA 

PARAM 

TPIMSH 

PARAM 


Description 

When MSTST=0, start iterations with a 
Phio distribution from a steady-state 
solution. When MSTST=0, U(I,J,K) is set to 
K and thickness effects are not included in 
the analysis. This is the "flat plate” 
analysis 

This variable along with ILAX determines 
the order of calculations; with ILAX=0 (row 
relaxation), ISWEEP has two possible 
values; (a)ISWEEP=0 indicates rows will 
be solved from the upper and lower 
boundary, alternating, in toward the wing; 
(b) with ISWEEP=0, rows will be solved 
from lower to upper boundary consecutively 
with ILAX=1 (column relaxation), ISWEEP 
has two possible values; (a) with 
ISWEEP=0, columns will be calculated 
starting from the trailing edge and moving 
forward (upstream) then coming back and 
calculating the nodes back of the trailing 
edge moving backward (downstream); 
(b)with ISWEEP=1 (with ILAX=1), 
columns will be calculated from upstream to 
downstream boundary 

With ILAX=0 (see also ISWEEP), 
relaxation using rows of points is used, the 
points forming a line parallel to the wing in 
an X-Z plane; with ILAX=1, relaxation 
using columns of points is used, the points 
forming a line perpendicular to the wing in 
an X-Z plane 

Constants required for convergence of 
row relaxation (See refs. 2 and 3.) 

Y coordinate of wingtip; scaled value 

Y coordinate of the inside edge of the 
control surface; scaled value 

This variable determines from what files 
the XYZ mesh size and the coordinates of 
each node are read. If TP1MSH=0, read 
IMAX, JMAX, KMAX, and the XYZ mesh 
coordinates from input. If TP1MSH=1, read 
them only from file named TAPEl. If 
TP1MSH=2, read them from both input file 
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Table 1. - (Continued) 


Card number Variable name Format 
2 to N (cont) 

LUPSYM PARAM 


IPLOT PRAM 

NDWING PARAM 

IMODE PARAM 


End of Namelist Variables 

Card number Variable name Format 

N+1 X 8F10.2 

to 

M 

M+1 Y 8F10.2 

to 

K 

K + 1 Z 8F10.2 

to 

L 


Description 

and file TAPEl. The data overwrites the 
TAPEl data. This last preserves file 
spacing on TAPEl 

If LUPSYM=0, perform a normal problem 
of full size. If LUPSYM=0 and the 
angle-of-attack THETA is 0, then assume 
the nodes above the wing are antisymmetric 
to those below and calculate only those 
below 

Not used 

The number of Y nodes on the wing 

If IMODE=l, a control surface mode is 
used, and the mode shape deflection is 
allowed only over the control surface 
defined by XA, XTE, and YSA. The axis of 
rotation is parallel to the Y axis at XA. For 
IMODE=2, a pitch mode is used, and the 
mode shape deflection is allowed over the 
entire wing. The axis of rotation is parallel 
to the Y axis at XA. When IMODE=3, 
flapping mode, mode shape deflection is 
allowed over the entire wing; the axis of 
rotation is parallel to the X axis at Y==AL 


Description 

The location coordinates of 

each of the IMAX X node locations; not 

input if TP1MSH=1; eight per card 

The location coordinates of each of the 
JMAX Y node locations; not input if 
TP1MSH= 1; eight per card (scaled 
values) ^ 

The location coordinates of each of the 
KMAX Z node locations; not input if 
TP1MSH=1; eight per card (scaled 
values) 
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Table 1. - (Concluded) 


Card number 

Variable name 

Format 

Description 

L+1 

XLE 

8F10.2 

The X location at the corresponding Y 

to 



node of the wing leading edge; NDWING 

P 



of these are input; eight per card 

P+1 

XA 

8F10.2 ' 

The X location at the corresponding Y 

to 



node of the aileron hinge; NDWING 

Q 



of these are input; eight per card 

Q+l 

XTE 

8F10.2 

The X location at the corresponding Y node 

to 



of the wing trailing edge; NDWING. 

R 



of these are input; eight per card 


5.8 PROGRAM OUTPUT 


5.8.1 PROGRAM RESULTS 

A listing of the output for a sample program is presented in appendix A. 

The printed output of the program consists of an initial block of information printing 
back the user’s input followed by information identifying the program options the user 
selected. 

The X, Y, and Z mesh data come next, read either from cards or binary file (TAPED. 
The mesh data are followed by the X locations on the wing of the leading edge, aileron 
pivot, or pitch axis and the trailing edge. Intermediate information regarding calculated 
variables and time used at routine calls will follow. This is followed by iteration prints 
giving data on how convergence to the required error difference is proceeding. 

If the program iterates to the maximum number of iterations specified (NMAX), the 
next data printed will be the complex Phii matrix. If the program stops before NMAX is 
reached, Phij will not be printed. If it is desired, a follow up run with NMAX=1 will 
always print the Phij matrix. 

Next, and finally, the pressure coefficients above and below the wing will be printed. 

Optionally, the user may print the pressure coefficients unextrapolated to the wing 
surface by setting variable NP to the iteration interval desired. 

5.8.2 PROGRAM DIAGNOSTICS 

"ERROR DIFFERENCE IS GREATER THAN 100 BETWEEN ITERATIONS. ERROR IS 
THE SOLUTION IS DIVERGING.” The program checks errors by calculatin g 
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the difference from iteration to iteration between corresponding nodes, saving the 
largest for comparison to the user-specified standard. The preceding message indicates 
the error is becoming larger too rapidly. 

"SOLUTION FAILED TO CONVERGE IN ITERATIONS lERR, 

JERR, KERR, ERRMAXl”. indicates the largest error found at the 

indicated XYZ node location is still larger than the user-specified standard and that the 
maximum number of iterations has been attained. 

"SOLUTION CONVERGED MAXIMUM ERROR IS ” indicates the user error 

standard has been reached and calculation will stop. 

"PLANE ERROR THE PLANE READ DOES NOT MATCH THE PLANE DESIRED J = 
— JT — lUNIT = — KAW= — ” indicates an error in the incoming Phij matrix on 
the indicated file. This error is probably caused by use of the wrong file for TAPE 10 
(Phij). 

"INFORMATIVE ERROR - - ISWEEP OPTION MUST BE 2 IF LUPSYM OPTION AND 
ROW OPTION SELECTED. ISWEEP SET TO 2” indicates user did not select correct 
option of ISWEEP for row relaxation. This also indicates ISWEEP was reset to the 
correct value. 

"STOP 1” - If this message appears in the dayfile, the program was unable to find the Z 
nodes just above and just below the wing. This may be caused by a Z node at Z = 0, 
which is prohibited. 


ORIGINAL PAGE IS 
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6.0 COMPUTER PROGRAM DESCRIPTION 


Flow diagrams for the program and its subroutines are listed in appendix B. 

6.1 OVERLAY STRUCTURE 

TEV174 consists of a (0,0) level overlay and four primary level overlays. The (0,0) 
overlay contains the program driver and several common usage subroutines. The (0,0) 
driving program called V174 does the following: 

1. Calls the (4,0) level overlay to input data and initialize variables and arrays 

2. Controls the number of iterations performed either by terminating because of 
convergence (or divergence) or terminating because the maximum number of 
iterations has been performed 

3. Calls the far-field boundary updating overlay (3,0) 

4. Calls the relaxation overlay specified by the user for row relaxation (1,0), 
specifying which X-Z plane is to be calculated 

5. Prints convergence data for each interation 

6. Saves the Phii matrix on file TAPElO after iterations have ceased 

7. Calls the final print package routine CPR for pressure differential printing 
The (4,0) overlay does the following: 

1. Reads inputs and prints them back 

2. Sets constants 

3. Calculates the body function 

4. Calculates the mesh dajta , , 

5. Initializes the steady-state Phio matrix with a constant or with data from another 
program 

6. Initializes the Phij matrix either with a complex constant or with data from a 
previous program run 

•••The (T,0) overlay updates an X-Z plane of the Phij matrix by column relaxation. It also 
causes the print display of the Phij matrix if the rhaximum number of iterations 
(NMAX) has been reached. 
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The (2,0) overlay, like the (1,0) overlay, updates the X-Z plane specified by the (0,0) 
overlay. It does this by row relaxation instead of column. It will also print the Phij 
matrix at iteration NMAX. 

The (3,0) overlay calculates and updates the nodes on the outside boundaries (farfield 
update overlay). 

6.2 COMMON BLOCK USAGE 

Table 2 describes usage of blocks of common variables. 

Boeing Commercial Airplane Company 
P.O. Box 3707 

Seattle, Washington 98124 
October 1977 


Table 2. ~ Blocks of Common Variables 


Block Usage 

CONST Control and calculation constants 

IMESH X, Y, and Z mesh constants 

WING Variables describing the wing 

CONST2 Control and calculation constants 

TIMES Print time variables 

CONE6 Row relaxation constant 

TROUB Troubleshooting variables 

ENTR Column relaxation control and calculation variables used in MATRXCF 

EPS Calculation constants 

FARE Far-field boundary calculation variables ■ 

IRELAX Relaxation variables used in MATRXCF and SRMATCF 

JAYS Phil incore plane designators. - Looking toward the wing from upstream; 

JL indicates the left X-Z plane. 

JC indicates the center X-Z plane. 

JR indicates the right X-Z plane. 

ITER Current iteration number and NMAX 

IVOL Wing integral variables 

MATRIX Matrices used to form the diagonals for solution; used in SRMATCF, 
MATRXCF, FOURDG, and TRIDIAG 

SUBER Contains two diagonals (see MATRIX) needed for PHICLM calculation 

used in SRMATCF, FOURDG, MATRXCF, and TRIDIAG 

SUBSUB Contains the fourth diagonal required for row relaxation; used in 
FOURDG and SRMATCF 


ORIGINAL 
OF POOR QUALITY 
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table.2. - (Concluded) . 

Block Usage 

PHI Contains the three X-Z planes currently being used by the program (see 

also JAYS) 

XYZ Contains the X, Y, and Z coordinates of the nodes at which Phij is to be 

calculated 

PRNT Not used, was a print variable 

IRLAXF Subsonic, supersonic node identification array 

ERRS Error convergence data 

10 File names 

FUL Upper and lower wing function 

DEL Phil array 

U Coefficient of <p term 

FXY Functions of the X, Y, and Z mesh 

BIAS Relaxation coefficients used to speed up convergence 

Note: The following common blocks are used only in the far-field boundary 
overlay (3,0). 

V Integral evaluation variables 

SUM Not being used in this version 
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.4363 

.41696-01 

.4363 

.41896-01 

.4363 

.41896-01 


.4363 

.41896-01 

.4363 

.41896-01 

,4363 

.41696-01 

.4363 

.41896-01 

.4363 

.41896-01 


"" • 436 3 

".41896-01 

• 4363 

■ ' .41696-01 







29 

.4363 

.44516-01 

• 4363 

,44516-01 

• 4363 

,44516-01 

.4363 

.44516-01 

• 4363 

.44516-01 




«A363 

.AA51E-0I 

.9363 

.9951E-CI 

.A363 

.AA51E-01 

.9363 

,995IE-C1 

30 .A363 

4 A 66 OE-OI 

.936 3 

.96606-01 

.A363 

,A660E-01 

.9363 

.96eOE-01 

.A363 

4 A 66 OE-OI 

.9363 

,9660E-01 

31 .A363 

.4817E-01 


.9817E-01 

.A363 

4A817E-01 

.9363 

, 9817E-C1 

.A363 

4A817E-01 

.9363 

.9617E-01 

32 4A363 

4A922E-01 

.9363 

.9922E-0I 

.A363 

.A922E-01 

.9363 

,9922E-01 

4A363 

,A922E-01 

.9363 

.9922E-CI 

33 .A363 

“T363‘ 

45027E-01 

.9363 

.5027E-C1 

450276-01 

.9363 

.5027E-C1 

4^36 3 

.5027E-01 

.9363 

.5027E-01 

3A 4«363 

451316-01 

.9363 

.5131E-01 

_ - - -,<>363 ' 

. .51316-01 

.9363 

.51316-01 

4 A363 

.5131E-01 

.9363 

.5131E-C1 

35 .A363 

.52366-01 

.9363 

. 5236E-CI 

*4363 

.52366-01 

• 4363 

.5236E-C1 

4^363 

.52366-01 

.9363 

,52366-01 

CP TIME TOT" 

22.796 SINCE 

LAST CALL 

.16C 

CP tine TOT* 

22.799 SINCE 

LAST CALL 

.0C3 

CP TIME TOT" 

22.996 SINCE 

LAST CALL 

.197 

N>IERRr JERt)>KERR> 

ERR MAXI 

1 39 ' 

I 3 

ITER 1 PLANE 

0 ERRMAXl" 

.60986-09 


ERRMAX • 49677E-0A ERROR ■ 

.10006-03 



SCLUTIDN CON'VERGtD.' MAXIMUM ERROR IS .60<»8£-0<i 
■~C> T'IM'E TOT" 25.04il SINCE LAST CALL Z.OVt 



■.9363'’ 

.■ 99516 - 01 " 

.9363 

.99516-01 

.9363 

,99516-01 

.9363 

.96606-01 

.9363 

.96606-01 

.9363 

.96606-01 

.9363 

.9660E-01 

.9363 

,9660 6-01 

.9363 

.96606-01 

• 9363 

.9617E-01 

.9363 

,9cl7E-01 

.9363 

.98176-01 

‘.9363 

.9617E-01 

.9363 

.9817E-01 

.9363 

.96176-01 

.9363 

.99226-01 

.9363 

.9922E-01 

.9363 

.99226-01 

.9363 

.99226-01 

.9363 

.9922E-01 

.9363 

.99226-01 

.9363 

.5C27E-C1 

.9363 

.5027E-01 

.9363 

.50276-01 

.9363 

.50276-01 

.9363 

.50276-01 

.9303 

.5027E-01 

.9363 

.51316-01 

.9363 

.5131E-01 

.9363 

.513a6-01 

.9363 

.51316-01 

.9363 

.5131E-01 

.9363 

.51316-01 

,9363 

.52366-01 

.9363 

.5236E-01 

.9363 

.52361-01 

.9363 

,52366-01 

.9363 

.52366-01 

.9363 

.5236E-01 


.60E-0A 



MACH NUHB6R • .86 

ckega'~« ,060“ 

amplitude • .026 

i ■ ■■■ ■ BELCH AlkFCIL 

PRESSURE COEFFICIENTS 
■ ABOVE AtRFCl'L 

A3QV6 - BELOW 


■ rr 

CHORD 


0.0000 X— 

leading edge* 

-i.cuoo 

AILERCN PIVOT. 

-1.0000 trailing 

ECG6< UOOOC 



10 

0. 

c. 

c . 

c. 

.2336*02 

-.5066*01 



11 


0. 

o"* 

c. 

.1606*02 

-.3936+01 



12 

0. 

c. 

0. 

c « 

.1506*02 

-.3236*01 



13 

0. 

c. 

0. 

c. 

,1396*02 

-.2836+01 

*” * * 

* “ * 

19 ‘ 

0. 

c * 

c. 

c. 

.1196*02 

-.2916+01 



15 

0. 

c. 


c. 

.1096*02 

• -.1966*01 



15 

0. 

c. 

0. 

c, 

.9276*01 

-.1606+01 



17 

6. 

c. 

c. 

c. 

.6916*01 

-.1916*01 



18 

0. 

‘c. - 

0, 

c. 

.6996*01 

-.1316*01 



19 

0. 

c. 

0. 

c. 

.8536*01 

-.102E+01 

. 

*** ” “ 

20 

0. 

c. 

^ • 

c. 

.8066*01 

-.6736*00 



21 

0. ' 

C, 

c. 

c. 

.8996*01 

-.5286*00 



Z2 

0« 

0* 

Oo 

0. 

.1276*02 

-.1586*00 



23 


C. 

0« 

C. - - 

,1236*02 

.2306*01 



29 

0. 

c. 

c. 

c. 

.6666*01 

.3676*01 



25 

0. 

0. 

Co. 

0. 

.3906*01 

.2578*01 


* 

26 

0. 

0. 

c. 

c. 

.2106*01 

.1958*01 



27 

0. 

c . 

c. 

c • 

. 1966*01 

.1578*01 



26 

0. 

c. 

0. 

c . 

.1066*01 

.1296*01 



29 

■ ' 0. 

c. 

'0. ■ 

c. 

.7656+00 

.1066*01 



30 

0. 

c. 

c. 

c. 

.6006*00 

.8836*00 



31 

0. 

c. 

0. 

c. 

.9696*00 

.7906*00 


— 

32 

0. 

0. 

0. 

c. 

.3826*00 

.6306*00 



33 

0. 

c. 

0. 

c* 

.2906*00 

.503E+00 



39 

0. 

c. 

0. 

c. 

,ie9£*00 

,3576*00 



"35“ • “ 

0. 

c. 

0. 

" c. 

0. 

0. 

AT 

CHORD 

Y( J )- 

.1701 X— 

LEADING edge- 

-loOOCO 

AILEFCN PIVCT* 

-1.0000 trailing 

6CGE> 1.0000 



io 

0. 

c. 

c* 

c. 

.2336*02 

-.5056*01 



11 • 

0. 

o’. 

Oo 

0. 

.1606*02 

-.3976+01 



12 

0* 

c. 

0. 

c. 

.1506*02 

-.322E+01 



13 

0. 


c. 

c. 

.1396+02 

-.282E+01 



19 

0. 

c. 

0. 

c * 

.1196*02 

-.2906+01 



15” " ■ 

0. 

c. 

' c. 

■ 0. 

.1096*02 

-.1956*01 



16 

0. 

c. 

0, 

c. 

.9256*01 

-.1596*01 



17 

0. 

c. 

u* 

c • 

.8896*01 

-.1906*01 



IT “ 

0. 

c. 

Oo 

0. 

,8936*01 

-.1306*01 



19 

0. 

c. 

0. 

0. ' 

.8526*01 

-.1016+01 



20 

0. 

0. 

0. 

c. 

,8076*01 

-.6626*00 

— — 


‘21 

0. 

c. 

c. 

c. 

.8966+01 

-.51 36+00 



22 

0. 

0. 

c. 

c» 

.1316*02 

-.8736-01 



23 

0. 

c. 

0. 

c» 

.1256+02 

•237E+01 



T9' 

■■’o; 

c. 

0. 

0. 

.6366*01 

.3656*01 



25 

0. ' 

0. 

0. 

0. 

.3166*01 

.2536*01 



26 

0. 

c. 

0» 

c. 

.1976+01 

.1926*01 

- “ 

— 

IT 

0. 

c. 

0. 

c. 

.1386*01 

.1556*01 



2B 

0. 

0, 

c. 

c« 

.102t*01 

.1286*01 


X 


-.960t«00 
-,9^0E+00 
-.9006*00 
-. 8 <. 0£*00 
-.7606*00 
-.6606*00 
-.6006*00 
-.5006*00 
-.600£'*00 
-.3006*00 
-.2006*00 
-.1006*00 
0 . 

.1006*00 

.2006*00 

.3066*00 

.*oce*co 

.5006*00 

.6006*00 

.7006*00 

.7606*00 

.6606*00 

.8606*00 

.9206*00 

.9606*00 

.1006*01 


-.9606*00 
-.9606*00 
-.9006*00 
-.8606*00 
T. 7606*00 
-.6606*00 
-.6006*00 
-.5006*00 
-.9006*00 
-.3006*00 
-.2006*00 
-.1006*00 
0 . 

.1006*00 

.2006*00 

.3006*00 

.9006+00 

.5006*00 

.6006*00 




OJ 

u> 


0 • 

Cm 

,757e+00 

.1056+01 

.7006+00 

0* 

Cm 

.5B06+00 

.8796+uO 

.7806+00 

c« 

c. 

.<>5BttOO 

.7336+00 

.8906+00 

G. 

c« 

.37Xe+00 

.6256+00 

.6806+00 

U* 

c. 

.202fc+OO 

.9996+00 

.9206+00 

C • 

c • 

.lB3e*00 

.3596+00 

.96 06 + 00 

0* 

Cm 

0. 

0. 

.1006+01 

l.COCO 

AILEPEh PIVCT- 

-1.0000 TRAUINC 

6D06* l.OCOO 


0* 

0 • 

.2326+02 

-.9996+01 

-.9806+00 

(r» 

c* 

.1596+02 

-.3936+01 

-.9906+00 

c. 

c. 

.1996+02 

-.3106+01 

-.9006+00 

c • 

Cm 

. 1336+C2 

-.2786+01 

-.8906+00 

G» 

Cm 

.1196+02 

-.2366+01 

-.7606+00 

0« 

Cm 

.1096+02 

-.1926+01 

-.6806+00 

0« 

Cm 

.9226+01 

-.1566+01 

-.6006+00 

0« 

c • 

.6866+01 

-.1376+01 

-.5006+00 

c* 

c* 

.8896+01 

-.1266+01 

-.9006+00 

0. 

c. 

.0986+01 

-.9756+00 

-.3006+00 

c. 

c* 

.6096+01 

-.6276+00 

-.2006+00 

c» 

Cm 

.8956+01 

-.9686+00 

-.1006+00 

c. 

Cm 

.1966+02 

.2556+00 

u. 

€• 

c. 

.1306+02 

.2696+01 

.1006+00 

c. 

c« 

.9776+01 

.3506+01 

.2006+00 

c • 

c. 

.2326+01 

.2366+01 

.3006+00 

0» 

c. 

.1576+01 

.1826+01 

.9006+00 

c • 

Cl 

.1176+01 

.1936+ul 

.5006*00 

c« 

c • 

.8966+00 

.1236+01 

.6006+00 

c* 

c« 

.6836+00 

.1026+01 

.7006+00 

* c* ’ 

c. 

.5306+00 

.6516+00 

.7806+00 

c« 

Qi 

.9106+00 

.7156+00 

.8906+00 

c » 

0* 

.3926+00 

.6106+00 

.S8C6+C0 

0. 

c« 

.2606+00 

.9886+00 

.9206+00 

c* 

€« 

.1706+00 

.3966+00 

.9606+00 

c. _ 

C. 

C. 

0. 

.1006+01 

l»GOC0 

AILCkCK PJVCT- 

'1.0000 TRAILING 

6DG6> ■ 1.0000 


0. 

C* 

,2316+02 

-.9696+01 

-.9806+00 

c« 

c* 

.1586+02 

-.3366+01 

-.9906*00 

Cf 

c* 

.1966+02 

-.3126+01 

-.9006*00 

G« 

Cm 

.1326+02 

-.'2726+01 

-.8906+60 

0« 

Cm 

.1186+02 

-.2316+01 

-.7606+00 

0« 

c. 

.1036+02 

-.1876+01 

-.6806+00 

Cm' 

c. 

.9166+01 

-.1516+01 

-.6006+00 

Cm 

c« 

.8806+01 

-.1316+01 

-.5006+00 

Cm 

c. 

.6826+01 

-.1206+01 

-.9006+00 

— c'.‘ 


■ .6906+01 

-..9126+00 

-.3006+00 

c« 

Cm 

.7986+01 

-.5696+00 

-.2006+00 

0« 

Cm 

.0906+01 

-.3936+00 

-.1006+00 

Cm 

Cm 

• ..976+02 

.9976+00 

0. 

c* 

C . 

,1266+02 

.2816+01 

.1006+00 

c • 

Cl 

.3986+01 

.3936+01 

_ .2006+00 

C » 

c. 

.1686+01 

"ii'iai+ol 

«300£+00 

0» 

G« 

.1306+01 

.1796+01 

.9006+00 

c* 

C 1 

.9986+00 

.1916+01 

.5006+00 

Cm 

Cl 

.7636+00 

.1166+01 

,6006+00 

Cm 

Cl 

.6076+00 

•9oXE>00 

.7006+00 



u> 

41 ^ 



30 

31 

32 

0. 

0. 

0. 

0. 

0. 

0_; 

c. 

0. 

c. 

0. 



sr 

39 

35 

i 

r CHORD Y(J». 

.6805 

X“ L6ADIK6 

EDGE 


lo 

0. 

0. 



11 

0. 

C. 



12 

.. 0. 

C. 



13 

6. 

- “c; 



19 

0. 

c . 



15 

0. 

c. 



16' 

0. 

c. 



17 

0. 

0. 



18 

0. 

c. 


i 

19 

0. 




20 

0. 

c. 



21 

0. 

0. 



22 

0. 

c. 



23 

0. 

c. 



29 

0. 

c. 



25 

0. 

■ o^. 



26 

0. 

c. 



27 

0. 

c « 



28 

0. 

c. 



29 

0. 

c. 



30 

0. 

0. 



Ji 

■ or. ■ 

o'. 



32 

0. 

c • 



33 

0. 

c. 




0^."‘ 

“c\ ■■ 



35 

0. 

c. 



T ChorO Y(JT^ 

■■.8506 

X— LEADING 

6DG6 


10 

0. 

C. 



11" " 

0. ■■ 

C. 



12 

c. 

C. 



. 13 

0. 

0. 



19 

0. 

"C. 



15 ' 

0. 

c. 



16 

0. 

c. 



ir' ■■■ 

-■ 0. 

c. 



18 

0. 

c. 



19 

0. 

c. 



20 . ■ 

0". ■ 

0. 



21 

0. 

c. 



22 

0. 

" ■ 

c. 



■23“ 


■ c. 



29 

0. 

c. 



25 

0. 

c. 



“2 ’6 

tr, - • 

cV ■' ■ 



27 

0. 

0. 



28 

0. 

c. 



29 ■ 

0. ■ ■ ■ 

c. 



30 

0. 

c. 








C. 

c. 

.9756+00 

. 8208+00 

0. 

c. 

.3776+00 

.6906+0C 

0. 

c. 

.3096+00 

.5906+00 

€• 

c. 

.2366+00 

.9726+00 

0. 

c » 

.1596+00 

.3356+00 

c. 


0. 

0. 

-1.0000 AI16RC6 PIVCT* 

-1.0000 TRAILING 

EDGE* l.OGOO 

c. 

c. 

.2296+02 

-.9756+01 

c. 

c. 

.1576+02 

-.3266+01 

t* 

c. 

.1976+02 

-.3026+01 

c. 

c. 

.1316+02 

-.2696+01 

0. 

c. 

.1176+02 

-.2236+01 

0. 

c. 

.1026+02 

-.1796+01 

0. • ■ 

c. 

.9066+01 

-.1996+01 

c. 

c. 

.8706+01 

-.1296+01 

c. 

G. 

.6706+01 

-.1126+01 

0* 

c. 

.8266+01 

-.8196+00 

0. 

c. 

.7876+01 

-.9826+00 

0. 

c. 

.8826+01 

-.2796+00 

0. 

c. 

■ .1956+02 

.8936+00 

0. 

c. 

.1176+02 

.2986+01 

c. 

c • 

.3006+01 

.3256+01 
■ .2l'3'6+’0i~ 

0. 

c. 

.1986+01 

0. 

0. 

.1076+01 

.1626+01 

c« 

c. 

.8536+00 

.1336+01 


1. 

.6856+00 

.1126+01 

0. 

L. 

.5386+00 

.9326+00 

0. 

c. 

.9296+00 

.7826+00 

0. 

c. 

.3386+00 

.6596+00 

0. 

c . 

.2776+00 

.5696+00 

0. 

c. 

.2126+00 

.9526+00 

o'. 

c. 

.1386+00 

.3216+00 

c. 

c. 

0. 

0. 

-iVOOOO" AILERCK PIVCT- 

-1.0000 TRAILING 

6CG6* '■ 1.0000 

0. 

c. 

.2266+02 

-.9576+01 

c. 

• ■ e. 

.1556+02 

-.3136+01 

c. 

c. 

,1956+02 

-.2906+01 

0. 

c. 

.1306+02 

-.2526+01 

0. 

0. 

.1156+02 

-.2l'2E+01 

0. 

c. 

.1006+02 

-.1696+01 

c. 

c. 

.6926+01 

-.1396+01 

G.' 

' C. ■ 

.859C+01 

-.1136+01 

c. 

0. 

.8516+01 

-.9926+00 

0. 

0. 

.8016+01 

-.6706+00 

0.'" 

c; 

.7696+01 

-.3326+00 

0. 

c. 

.8716+01 

-.8796-01 

c. 

0. 

.1526+02 

.1756+01 

_ .. . 

■ c. ‘ 

.1076+02 

' .3306+01 

0. 

c . 

.6726+00 

.2516+01 

c. 

c. 

.7136+00 

■ .1756+01 

0« 

"c. 

.7716+00 

.1916+01 

c. 

c. 

.6996+00 

.1206+01 

c. 

c. 

.591E+0C 

.1036+01 

— "o;' ■ 

r.’ 

,9756+00 

.3696+00 

0* 

c. 

,3786+00 

.7396+00 


.760E+00 

.8^0fc+C0 

.6806^00 

. 9208^00 

.960E«00 

.10bt>01 


-.9806+00 
-.9906+00 
-.90Ct+00 
-.8906+00 
-.760E+00 
-.6806+00 
-.6006+00 
-.6006+00 
-.9006+00 
-.3006+00 
-.2006+00 
-.1006+00 
0 • 

.1006+00 
.2006+00 
' ;‘3006+00 
.9006+00 
.3006+00 
.6006+00 
.7006+00 
.7806+00 
.8906+00 
.8606+00 
.9206+00 
.9606+00 
.1006+01 


-.9806+00 
-.9906+00 
-.9006+00 
-.8906+00 
■■ -.’7606 + 00 
-.6806+00 
-.6006+00 
-.3006+00 
-.9006+00 
-.3006+00 
■■-.2006+00 
-.1006+00 
0 . 

■ .1006+00 
.2006+00 
.3006+00 
.9006+00 
.5006+00 
.6006+00 
.7006+00 
.7806+00 




Bp 

I? 


31 

0« 

c. 

3Z 

0, 

c* 

33 

0. • 

c. 

- 3^- 

0. 

C. 

35 

0. 

c. 

“AT CHORD yO'i- 

1.0207 X — 

LEADING EDGE' 

10 

0. 

C. 

'^rx' 

0. 

C. 

IZ 

0* 

G. 

13 

0. 

c. 

14 

0* 

c* 

15 

0. 

C. 

16 

0. 

c. 

17 ■ 

0. 

0. 

le 

0» 

c » 

19 

0. 

G. 

20 

Oi 

G. 

21 

0. 

G. 

22 

0. 

C. 

23‘ 

Oo 

0. 

24 

0. 

c. 

2^ 

0* 

t* 

■ “ 26 

0. 

c* 

27 

Oi 

c. 

2B 

0. 

0. 

- 29 

0* 

G* 

30 

o« 

c. 

31 

0* 

c. 

32 

0« 

C. 

33 

0. 

c. 

S'. 

0. 

c. 

35' 

0. 

c. 

A1 CHORD 7(J)- 

1.190B X— 

■ LEADING edge 

10 

0. 

c. 

11 

0* 

0. 

^2 

0* 

€• 

13 

0. 

C. 

l<i 

0. 

c. 

15 


0« 

16 

0. 

c« 

17 

0. 

c. 
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c* 
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c. 
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c. 
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c« 
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.3026+00 

.6216+00 

.89CE+00 

1 2966+00 

.532t+0C 

.atOE+GG 

.1906+00 

.9276+00 

.9206+C0 

.I23b+0C 

• 30ti6+00 
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0. 
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-.932t+01 

-.96CE+C0 

.1526+02 
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.1936+02 

-.2796+01 

-.9006+00 

.1276+02 

-.2376+01 

-♦690t+00 

.1136+02 
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-,76CE+0Q 

.9606+01 

-.1566+01 
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.6696+01 
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.8316+01 
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-.6226+00 

-.9006+C0 
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-.1296+00 

-.2006+00 

.6636+01 

.190E+0C 

-.lOOE+CO 

.1976+02 

.3136+01 

0. 

.6556+01 

.3956+01 

.iOOt+CO 

-.1066+01 

.1136+01 

.2G06+G0 

.9516+00 

.130E+G1 

.3006+C0 

.6776+00 

.1206+01 

.9006+CG 

.6386+00 

.1076+01 

.5006+00 

.5936+00 

.9916+00 

.6006+00 
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.6056+00 

.700E+CO 

.3976+00 

.6896+00 
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-.2796+01 

-.9906+00 

.1396+02 

-.2526+01 

-.9006+00 

.1236+02 

-.2176+01 
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-.6806+00 
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-.1056+01 
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.6766+01 
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-.3006+00 

,6976+01 
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0. 
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.3006+00 

.6096+00 

.1126+01 

.9006+00 

"■ " .6636+00 

,991^+OC 

.5006+00 

.5336+00 

.8706+00 

.6006+00 

.9176+00 

.7966+00 

.7006+00 

.3266+00 

.6356+OG 

.7806+00 

.2576+00 

.590E+OC 

.8906+C0 
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c. 
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30 
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c. 
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c. 

c. 
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c. 
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0. 

c. 

0. 

0. 
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c. 

c. 

.297E«02 

-.361E+01 

-.9806+00 

0. 

c. 

.lA2ei02 

-.2Aae+01 

-.9A0E+0C 

c. 

c. 
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-.225E+01 

-.9006+00 

c. 

c. 
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-.191E+01 

-.8A0E+00 

c . 

c. 

.103E>02 

-.15AE+01 

-.7o06+00 

0. 

c. 

.079E+O1 

-.1156+01 

-.6806+00 

0. 

c * 

.765E+01 

-.E20E+00 

-.6006+00 

L. 

c . 

.703EA01 

-.572E+00 
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c. 
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.652E+01 
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-.A006+00 

u. 

c. 

.567E+01 
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c. 

c. 
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lI 
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-.1006+00 

0. 

c. 
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0 . 

c. 

c. 
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c. 
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0. 

c. 
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c. 
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.1026+01 
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.2976+00 

.5796+00 

.7806+00 

c. 

c. 

.232E+00 

.A93E+00 

.dAOE+CO 
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c. 
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.8806+00 
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c. 

c. 
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c. 
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.1926+02 

-.3076+01 

-.9806+00 
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• 131E+02 

-.2086+01 

-.9A06+00 
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.1226+02 

-.1896+01 

-.9006+00 
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. 1066 + 02 

-.1566+01 

-.8A06+00 

c. 


.9006+01 

-.1196+01 

-.7606+00 

c. 


.7A96+01 

-.6276+00 

-.6806+00 

■ c'.~ 


.6A7I+01 

-.5536+OC 

-.6006+00 

0. 


.5636+01 

-.2786+00 

-.5006+00 

c. 


.+656+01 

-.223E-01 

-.9006+00 

c. 


.9366+01 

.1896+00 

-.3006+00 

c. 


.A676+01 

.A186+00 

-.2006+00 

c . 


.568E+01 
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-.1006*00 

C. ■ 


.9076+01 

.1516+01 

0. 
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.lAAE+01 

.1206+01 

.1006+00 
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.1196+01 

.1066+01 

.2006+00 
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.92AE+00 

.9A8E+00 

.3006+06 

G. 


.7216+00 

.6536+00 

.9006+00 

c. 


.5616+00 
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c. 
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33 
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» 79e>£ +01 

-*101£+tl 
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c* 
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c • 

c* 
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c • 
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c . 
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c* 
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c . 
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c* 

c« 
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0 • 

G* 
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G * 
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c* 
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c < 
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c* 

c* 
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^♦f»CiC t ♦Cu 
-.400k+CO 
“* 3uGk^C0 
-•^OOfctCO 
-.lOCt^CO 

c« 

t lCCt'^00 
»2uOb*CC 
, 300[400 
*<»UV;b«UU 
•iUCt+CO 

• bOCrOCC 
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, d 0 C t ♦ C 0 
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APPENDIX B 








PROGRAM/SUBROUTINE 
FLOW DIAGRAMS 


PROGROH TEV17U aNPUT=1002.0UTPUT»l002,TfiPE7,DEBUG.TfiP£5*lNPUT, 
R TRPES-OUTPUT. TflPEl ■= 1 002, TRPEl 0- 1002, 
7fiPElUJ0O2.TflPEI2=1002.TflPEI3»l00a 


PVIRPflSE 


0VERRLL PROGfWH PURPOSE 

TEV1714 INCLUDING OVERLRTS COMPUTES THE UNSTERDT TRRNSONIC 
TRRNSONIC PERTURBATION FLOH ABOUT fl KARMONICALLT 
OSCILLATING THREE DIMENSIONAL HING. USING FINITE DIFFERENCE 
TECHNIQUES. 

OVERLAT PURPOSE 


OVERLAY 0.0 DOES THE FOLLOHINGs 


SERVES AS THE OVERALL DRIVER, , 

SERVES AS THE INITIRLLIZHTION VARIABLE STORAGE AREA, 
CONTROLS THE ITERATION COUNTER AND TERMINATES AS NECESSARY, 
CHECKS THE ERROR DIFFERENCE BETWEEN ITERATION TO SEE IF„ 
CONVERGENCE OR DIVERGENCE HRS OCCURRED AND TERHINATES IF 
KECESSflRY 

CONTROLS THE FRRFIELB BOUNDARY UPDATING, 

CONTAINS UTILITY ROUTINES NEEDED BT THE OTHER OVERLATS. 


PROGRAM ROUTINE TEV17U PURPOSE 

THE TEV174 ROUTINE IS THE OVERALL DRIVER AND PERFORHS R TO E RBO 


AUTHOR T.D BUTLER 
REFERNCES 


LANGUAGE FTN 4.S 420 


1. EHLERS.F. EDWARD. A FINITE DIFFERENCE METHOD FOR THE SOLUTION 
OF THE TRRNSONIC FLOW AROUND HARMONICALLY OSCILLATING 
WINGS. 

(BOEING NASA CONTRACT .NASA CR-2257) 

2. WERTHERILL, WARREN H.. EHLERS.F. EDWARD. SEBASTIAN. JAMES C.. 
COMPUTATION OF THE TRANSONIC PERTURBATION FLOW FIELDS 
AROUND TWO AND THREE DIMENSIONAL OSCILLATING WINGS. 
(BOEING-NASR CONTRACT. NASH CR“2599) 

DESCRIPTION 

PROGRAM ROUTINE TEVnH HILL pErFORM ItE TASKS IN THE FOLLOWING 
ORDER: 

A] CALL THE INITIALIZATION OVERLAY FOR INPUT AND OTHER INITIALIZ 
B1 START THE ITERATION LOOP. IT HILL LOOP NMAX TIMES MAXIMUM, 

Cl IF INDICATED FOR THIS ITERATION CALL THE FRRFIELD UPDATE 
OVERLAY. 

D) DEPENDENT ON INPUT VARIABLE ILAX (SEE PflOGRAH INIT) 

CALL THE ROW RELAXATION OH THE COLUMN RB.AXTION CALC- 
ULATION OVERLAY (SEE REFERENCE 2) 

El WRITE ON TAPE THE FINAL XZ PLANE FOR THIS ITERATION, 


Figure B-1 .—Program VI 74 
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c 

c 

c 

c 

c 

c 

c 

c 

c 
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n TEST Fca CONVERGENCE OR DIVERGENCE. TERMINATE THE ITERATION 
LOOP IF EITHER IS, PRESENT, 

END OF THE ITERATION LOOP ' 

GJ SAVE THE PHI I DISTRIBUTION ON FILE TAPEIO, 

HI OUTPUT THE FINAL PRINTS 
I) EXIT 

NOTE SUBJECT TO CONVERGENCE/DIVERGENCE CRITERIA DO C TO F 
NMAX TIMES. 

MODIFICATION HSITORT. 


TEVITN 3-D TRANSONIC SHALL DISTURBANCE FLOW OVER WING 
VERSION 3-13-75 THREE 0 OUT OF CORE SOL 
VERSION 3-15-76 LO-UP STHMETRT-TB 
VERSION 5-3-76 THREE-D ROW RELAXATION-TB 
UNSTEADY TRANSONIC SMALL DISTURBANCE FLOW OVER AIRFOILS 

INPUT 


INPUT OF THE INITIAL BCD AND BINARY DATA IS ACCOMPLISHED IN 
THE INITIALIZATION OVERLAY t PROGRAM INIT IIIDL.il.O) 

OUTPUT 


THIS ROUTINE OUTPUTS THE SAVEABLE AND REUSABLE PHIl 

t 

DISTRIBUTION IN BINARY FORM ON TAPEIO. ' 

ALSO IT OUTPUTS PRINTED ERROR MESSAGES AND CONVERGENCE 
DIVERGENCE MESSAGES. 

OTHER I/O IS ACCOMPLISHED BY THE OVERLAYS AND ROUTINE 
CPfl IN THIS OERLAY. 

THE OTHER OVERLAYS MAY CALL THE UTILITY ROUTINES IN THIS 
OVERLAY FOR I/O. 


FSMACH 

MACH NO. 


DELTA 

THICKNESS RATIO 


THETA 

ANGLE OF ATACK 


OMEGA 

ANGULAR REDUCED FREQUENCY 


GAMMA 

RATIO OF SPECIFIC HERTS FOR AIR 

AL 

HINGE LINE 


IHAX 

MAX X NODE COUNT 


UMAX 

MAX T NODE COUNT . 


KMAX 

MAX Z NODE COUNT 


IS 

STARTING X NODE LIMIT FOR VOL 

INTEGRAL FOR MING 

IE 

ENDING X NODE LIMIT FOR VOL INTEGRAL FOR WING 

KS 

STARTING T NODE LIMIT FOR VOL 

INTEGRAL FOR MING 


Figure B-1.— (Continued) 
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c 

c 

c 

c 

c 

c 

c 

c 


KE ENDING T NGDE LIMIT FOR VOL INTEGRAL FOB KING 

NHflX HAXIHUH POSSIBLE COUNT OF ITEHRTI0N3 

Nfl FRR FIELD UPDATE CTCLE COUNTER. 

UPDRTE FRRFLD AFTER EACH NR ITERATIONS 
ERROR ERROR TARGET. WHEN ERROR IS LESS THAN THIS STOP 

NP PRINT ITERATION COUNT 

INC START RESTART VARIABLE. INC^O ON START. 

ORF SUBSONIC BIAS TO ACCELERATE CONVERGENCE,^ 

URF SUPERSDNC BIAS TO ACCCELERATE CONVERGENT 

HSTST when equal TO ZERO START WITH DBATB FROM 

ISHEEP DEFINES CRDER^oPnOBE^SOLUTION. PARTICULARLY ROWS. 

NVOL counter for VOL INTEGRAL INaUSION IN FRRFLD 

UPDATES 

ILRX RELAXATION PROCEDURE SELECTOR. NOT USED PRESENTLY 

CONPXT constant for PHI XT 

CONES constant REQUIRED FOR CONVERGENCE OF RON 
RELAXATION 

YT Y COORDINATE OF WING TIP 

YSA T COORDINATE OF INSIDE EDGE OF CONTROL SURFACE 


PERFORM INPUT AND INITIALIZATION 


I 


CALL OVERLAY (SLIIlDL.tl.Ol 


Q 


I 


THIS IS THE RELAXATION ITERATION LOOP 




I 


NFLIP=0 

I 


DO SCO ITER=l,NMflX 

T 



7879 CONTINUE 

CALL SECOND IT I HE) 
ETIME=TIHE-TIHEP 
TIMEP=TIME I 

KRITE (6.9100) TIME.ETIHE 
PRINT 7878 . N 






Figure B~1.— (Continued) 

























Figure B-1. -(Continued) 














Figure B-1.— (Concluded) 
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\ SUBROUTINE H1SP^(IH. JM.KMIX.N.JI 7 


PURPOSE 

ThIs'rOUTINE displays one plane of the COHPlfX mi 

ONE DISTRIBUTION AT EACH CALL. THAT IS ONE J PLANE 
tXZ PLANE AT A T NODE POSITION ON THE f AXIS! 

Z IS INCREASING BOHN THE PAGE. WHILE X IS INCREASING 
LEFT TO RIGHT. USUALLY Y IS INCREASING FROM PLANE TO 
^NE BUT W IS BEPENDENT ON THE CALLING ROUTINE 
CURRENTLY THE CALLING ROUTINES ARE SET FOR AN INCREASING T. 

AUTHOR T.D. BUTLER G25N11 HS/TS-SI PH 237957D 
LRNGUHGE FTNN.6 N20 

INPUT 

VARIABLE HETHQD BESCRIPTION 


IH 

JM 

KM 

N 

J 

PHI 

X 


CALL PAR MAX COUNT OF I THE X NODE COUNTS 

CALL PAR MAX COUNT OF J THE Y NODE COUNTER 

CALL PAR MAX COUNT OF K THE Z NODE COUNTER 

at }BI 5SK»igTF'?SE"SBNT J PL»NE 

PUM PLUS THE L»ST 

COHH/JATS/™C iPtHE NODE NUMBER IN THE PHYSICAL 


I 


PHI MATRIX IN CORE OF THE CURRENT XH PL 
JL AND, JR ALSO PRESENT BUT NOT ^ 

USED ARE THE NODE NUMBERS IN THE IN COR 
PHI MATRIX OF THE RESPECTIVELY. 

LAST AND NaT PHI PLANES. 

OUTPUT 

THirROUTINE OUTPUTS. AS DESCRIBED IN THE PURPOSE ABOVE. 

THE JTH XZ PLANE OF THE PHI 1 DISTRIBUTION IN BCD 
PRINT FORM. 



200 CONTINUE 
HR ITE (6.2021 
WRITE (6. 102) J 



ORIGINAL 
OF POOR 


WRITE (6,202) 
JPM=J+I 
WRITE (6, 1021 JPH | 

, t: i" 

DO U50 IslI.IM 

1 


WRITE 16,2011 1. (PHI (1. JR.Kl ,K=t,KMlXI 


Figure B-2.Subroutine DISP 







|H50 CBNTINUtI 


500 CONTINUE 
RETURN 


Figure B-2.— (Concluded) 





\ SUBROUTINE PLqNEtJJ.JL.JC.JFIl X 
PURPOSE 


THIS ROUTINE INPUTS ONE PLANE OF THE PHI 1 MATRIX AT A 
TIRE ADJUSTING THE JAYS COMMON BLOCK TO INDICATE HHim 
OF THE THREE IN-CORE PLRNES IS LEFT. CENTER. AND RIGHT. 

fflJTHOR T.D. BUTLER GaSUN MS 73-3>l 
LANGUAGE FTN N.6 U20 ’ 

INPUT 

VARIABLE METHOD DESCRIPTION- 

.. CALCULATED 

JL ‘ CALL PAR LEFT PLANE IN CORE CO^ TO JJ-1 
JC CALL PAR CENTER PLANE OF THE INCOHE PLANE 

CORRESPONDS TO THE JJ PLANE - 
JR Cffl-L PAR RIGHT PLANE NODE OF INCORE PHIIS 

CORRESPONDS TO JJ+l 





Figure B-3.Subroutine PLANE 




















Figure B-3.~(Continued) 










Figure B-3.— (Concluded) 


50 



ork^al page is 

OF BOOK QUAliW 



\ SUBROUTINE CPR tlPLOTJ / 

1 


PURPOSE 

THIS ROUTINE MODIFIES FIND OUTPUTS THE PRESSURE COEFFICIENTS 


(Wl=TH£Tfl/57.29S78 
WRITE (6,90001 FSHRCH.OMEGfl.RMl 
EPS2=-2 , kEPSkS 7. 29578/THETfl 
REWIND 11 
REWIND 12 
REWIND 13 
COM=IMRGNOHEGfl 
X=1 


I 


1300 CONTINUE 

I 


DO 200 J=l,JSl l 
^ 



Figure B-4.^Subroutine CPR 
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X SUBROUTINE mi U. JC] / 


PVJRPQSE 

CfiLCULfiTES THE PRESSURE fiCROSS THE AIRFOIL 
PRESSURE JUMP RCROSS AIRFOIL ' 




Figure B-5.— Subroutine DELT 
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Figure B-6.— Complex Function DELPHI 
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■ Figure B-7.— Program IN IT 
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Figure B-7.— (Continued) 
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Figure B-7.— (Concluded) 














Figure B-8.— (Continued) 
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Figure B-9.— Subroutine MESH 
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Figure B-W.^Subroutine INCOND 
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Figure B- 12.— Program COL LAX 





















Figure B-1 2.— (Concluded) 
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Figure B-73.— Subroutine TRIDIAG 
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Figure B-14.^Subroutine MATRXCF 
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Figure B- 14. —( Continued) 





Figure B- 14.— (Concluded) 
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Figure B-15.— Program ROW LAX 
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Figure B-1 5.— (Continued) 




Figure B-1 5. —(Concluded) 
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Figure B- 16 .— Subroutine FOURDG 
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Figure B- 17.— Subroutine SRMA TCF 











Figure B-1 7.— (Continued) 
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Figure B-1 7.— (Continued) 
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Figure B-7 7.— (Continued) 
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Figure B-7 7.— (Concluded) 
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Figure B-18.— Program FARFLD 
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Figure B-18.— (Continued) 
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Figure B-1 8.— (Continued) 







Figure B-1 8.— ( Continued) 



Figure B-1 8.— (Continued) 
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Figure B-18.-(Continued) 


















Figure B-1 9.— Complex Function FCNX 
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